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The Kirchhoff-Braess Paradox and Its Implications for Smart Microgrids 

John Baillieul, Bowen Zhang, Shuai Wang 

Abstract 

Well known in the theory of network flows, Braess paradox states that in a congested network, it may happen that adding 
a new path between destinations can increase the level of congestion. In transportation networks the phenomenon results from 
the decisions of network participants who selfishly seek to optimize their own performance metrics. In an electric power 
distribution network, an analogous increase in congestion can arise as a consequence Kirchhoff’s laws. Even for the simplest 
linear network of resistors and voltage sources, the sudden appearance of congestion due to an additional conductive line is a 
nonlinear phenomenon that results in a discontinuous change in the network state. It is argued that the phenomenon can occur 
in almost any grid in which they are loops, and with the increasing penetration of small-scale distributed generation, it suggests 
challenges ahead in the operation of microgrids. 


Index Terms 

distribution networks, network congestion, loss cost, price of anarchy 


I. Introduction 

A good deal of current research on the operation of smart grids has been focused on the information structures and protocols 
that enable operation, [8],[9],[10],[11],[12]. To facilitate load management at the level of a residential or commercial building 
microgrid, the authors have proposed the concept of packetized direct load control (PDLC). Here the term packetized refers 
to a fixed time energy usage authorization, with the emphasis being on usage by thermostatic appliances (water heaters, 
refrigerators, air conditioners, etc.) that have been aggregated into pools defined by common operating characteristics. Taking 
HVAC, for example, consumers in each room in a building choose their preferred set point, and then an operator {smart 
building operator, SBO) of the local appliance pool will determine an appropriate comfort band around the set point. Our 
previous work has shown that PDLC is capable of ensuring consumer comfort while at the same time reducing power 
oscillations that occur when no control is applied and appliances operate independently, [8] . 

While packet-switched energy distribution protocols appear to have numerous advantages in controlling demand oscillations 
from predictable loads, open questions remain concerning the operation of small-scale distribution networks in which multiple 
distributed generation (DG) sources inject power with hard-to-predict intermittency. In such settings, stable and secure grid 
operation will be increasingly challenged by non-radial distribution network topologies in which lines connecting loads and 
power sources are opened and closed according to arrays of factors including weather, time-of-day, and real-time energy 
market conditions. 

Focusing the discussion on network topology control. Section II begins with a review of Braess paradox in transportation 
networks. Section III focuses on an analogous phenomenon in a voltage controlled circuit. Here it is shown that changing 
the circuit topology by adding a small load can lead to relatively large losses in the circuit as a whole. Section IV discusses 
similar sensitivity in optimal power flows within a simple non-radial distribution grid. Section V extends the analysis of 
Section III to voltage-controlled circuits of arbitrary size. Open problems are discussed in Section VI. 

II. Problem Formulation 

The principal challenge addressed here is the management of network congestion in the presence of load and DG 
(distributed generation) uncertainties. We begin by recalling the well known Braess paradox that is generally associated 
with congestion in transportation networks. The setup, in simple form, is shown in Fig.[T] There is a network (of roads) with 
an origin O and destination D. A certain number of travelers will make the journey, and in the network in Fig. [ijA) they 
have a choice of the route with segments A® or segments CD. Congestion may enter either route in terms of travel times 
that depend on the number of users traveling on each segment. If / denotes the number of voyagers on the segment, the 
travel times on segments A and D are the same and equal to / + ,8 for some constant /3. The travel times on segments ® and 
C are similarly equal to af for some positive a. There are many different values of the parameters used in the literature, 
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but the basic idea is that because the left-hand and right-hand routes in Fig. [^have the same congestion cost, (a -b 1)/ + /3, 
introducing the cross link will break the cost symmetry and could cause the cost of travel to increase. Taking the particular 
values of [6]: a = 10,/3 = 50, 7 (-) = 0, and letting the total number of travelers be / = 6, we find that without the cross 
link, the best choice for minimizing travel time is for three of the travelers to choose the left-hand route and three to choose 
the right hand route. This is a Nash equilibrium. The travel time for each traveler is {a + l)f + f3 = 11 • 3 -|- 50 = 83. When 
the no-cost cross link is present,however, travelers will observe a possibly shorter route given by following the segments 
F-cross-link-!B. Indeed, if only three of the six travelers took this route, the travel cost could be as low as 2af = 60. 
Unfortunately, all six travelers may choose the route, in which case, the cost becomes 120. Taking the cost to be travel-time, 
Braess paradox is that adding a delay-free travel link can actually increase congestion and the users’ travel time. 




(A) (B) 


Fig. 1. The classic Braess paradox of congested network flows. 

Braess’ paradox in this setting involves selhsh social choices. There is a similar apparent paradox in electric circuits, as 
noted by Cohen and Horowitz ([4]) that is a consequence of the laws of electrophysics. In Fig. the question immediately 
arises as to how the horizontal connection changes the characteristics of the circuit on the top and bottom links. Specihcally, 
suppose that R 2 > Ri- Elementary invocation of Kirchhoff’s current and voltage laws indicates that in the absence of the 
horizontal link, the currents ii and 12 will be equal, but if the link is added with a moderate value of the resistance R 3 , 
we will have ii > 12 - This change is not surprising and is consistent with the observation in [4] that introducing the path 
changes the voltage drop across the circuit. It is also consistent with the observations in [2],[3] that adding this link may 
worsen congestion in similar models of power grid interconnections. The interesting question posed in [2] and [3] is whether 
(and when) there is a useful tradeoff in grid design that balances the increased reliability of larger numbers of power lines 
against the congestion that may occur due to adding these lines to the network. In what follows, we study a related question 
of when congestion may unexpectedly occur in a distribution network under the PDLC protocols described in [8]-[12]. 

Historically, distribution networks have had radial topologies, with distribution lines leading from a single trunk or power 
source to various commercial and residential loads. This standard topology is likely to change with the ever increasing 
penetration of distributed generation in the form of wind, solar, plug-in electric vehicles, and many other forms of alternative 
energy. Microgrids as depicted in Fig. will become prevalent, and these will increasingly resemble miniaturized transmission 
networks in which sources and loads will be connected through a multiplicity of line links that can be opened and closed 
as needed to maintain the needed balance between capacity and demand. It is against this backdrop that we examine the 
question of how the increased penetration of spatially distributed generation together with new electricity market models 
aimed at managing demand response will challenge the operating security of microgrid distribution networks. 

HI. The effect of DG intermittency on congestion sensitivity to small changes in load 

Recent work has demonstrated how various communication protocols can be effectively employed in networks of smart 
microgrids to improve energy efficiency, decrease demand volatility, and ensure customer satisfaction—[8],[9],[10],[11],[12]. 
Although this work has examined packet switched energy delivery in terms of temporal uncertainty in the operation of 
microgrids, the effects of spatial and network topology variability and uncertainty of both demand and DG capacity is not 
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Fig. 2. Distribution grids of the future will feature increased amounts of local generation involving power sources with widely varying capacities and 
increased operating uncertainties. (Figure from Arghandeh et ai [1]). 


well understood. While the precise magnitude of costs associated with mitigating the uncertainties of renewable generation 
sources is not known, some estimates suggest that higher reserve margins will be required. For instance, the historical 
averages of reserve requirements in the power grid point to 7% to 8 % as generally sufficient to handle contingencies. There 
are now predictions that if renewable penetration gets to the 33% level (still a long way off) these requirements may go as 
high as 15%, [7]. 


Ei 

E2 



Fig. 3. A simple circuit with DC-voltage sources and resistive loads. 

The implications of renewables for transmission networks remains a work in progress. With the increasing likelihood that 
distribution networks will incorporate small-scale distributed generation as depicted in Fig. we turn to the question of 
how these intermittent sources will affect network congestion on a small scale. Our working definition of congestion will 
be in terms of Fig. where we say that there is congestion if there is a significant difference between cutTents ii and i 2 
and more importantly, between the energy losses due to heat ifRi and We note that if there is no horizontal link (i.e. 

if i 3 = 0 or equivalently, if i ?3 ^ c»), then ii = 12 , and if i?i = R 2 , the heat losses at each resistor are equal as well. 
These conclusions remain true irrespective of the magnitudes of the voltage sources. If i?i = R 2 and the voltage sources 
happen to be equal in magnitude, then ^3 = 0 no matter what value is assigned to R^. If there is any imbalance in the 
voltages El and E 2 , or if one of them is zero (think of a wind turbine or solar array being out of service due to weather 
conditions), then a small value of R 3 on the cross link of Fig. can produce a very large difference in the cutTents ii and 




























4 


12 - We summarize this in the following: 

Proposition 1: Suppose Ri = R 2 = R and that Ei = 0. Then if the cross link is not connected in Fig. we have that 
*1 = *2 = E 2 / 2 R while if the cross is connected and R 3 is positive but small, we have 0 ~ ^ t 2 E 2 /R. Moreover, 

the total heat loss across across the entire circuit is approximately twice the loss if the cross link is disconnected. 

Proof: A simple application of Kirchhoff’s circuit laws shows that the currents ij are related to the circuit voltages 
and resistances by relationships that are linear in the voltages but subtly nonlinear in the resistances. 


ii 

*2 

*3 


/ El (R2+R3)+E2Ft3 \ 

Ri R2+Ri R3+R2R3 


Ri R2+Ri R3+R2R3 


—E1R2+E2R1 

Ri R2+Ri R3+R2R3 


/ 


( 1 ) 


The proof follows by setting Ri = R 2 = R and considering the equation for small positive values of R 3 . ■ 

The dramatic effect on current flow produced by connecting the cross-link in Fig. is illustrated in Fig. Connecting 

the link with a very small resistance i ?3 brings the current ii to essentially zero, while doubling the current 12 - It is not 
surprising that the losses on the line through which 12 flows increase dramatically. It is also the case that the total losses 
for the circuit double when the connection is made. 



Fig. 4. The effect of connecting the cross link in Fig.|^ The magenta line indicates the per unit current values * 1,12 when the cross is not present. The 
red, blue, and green curves are the values of the currents ii,j 2 , and 13 (as functions of the resistance iJa) when the cross link is connected. R\ = R 2 = R 
in all cases. 


Definition 1: The Loss Cost of the Link (LCL) is defined as the ratio between the total network losses before and after the 
addition of a new link or other capacity enhancement. Denoting the system losses before and after the capacity enhancement 
as Loss' and Loss, respectively. The LCL is defined as 


LCL = 


Loss 

Losf 


( 2 ) 


The LCL concept is our version of the price of anarchy that is discussed in routing congestion problems in transportation 
networks, [ 6 ]. Proposition 1 shows that the LCL for the balanced load network with Ri = R 2 is 2. The following Corollary 
extends the result to imbalanced networks of parallel connected resistors. 

Corollary 1: Suppose Ei = 0, if the cross link is connected and R^ is positive but small, then the price of anarchy of 
the Kirchhoff-Braess network is R 1 /R 2 + 1. 
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Proof: Based on Kirchhoff’s circuit Law, the current ij that flows through the resistor Rj,j = 1,2,3 are 


*1 = 

*2 = 

*3 = 


E 2 R 3 


R1R2 + R2R3 + R3R1 
E 2 {Ri + R 3 ) 

R 1 R 2 + R 2 R 3 + R 3 R 1 ' 
E 2 R 1 


R 1 R 2 + R 2 R 3 + R 3 R 1 

The total system loss for small value of i ?3 i?i, R 2 is 


Loss = i'jRj 
i=i 


El 

R 2 


On the other hand, the system loss before adding R 3 is 

Loss' = 

From the above two equations, we have the LCL is 


El 


(i?l + R2) 


(3) 


(4) 


LCL 


Loss 

Loss' 




(5) 


This ends the proof of Proposition 1. ■ 

A general scenario for Proposition 1 is to consider that two voltage sources are connected in the opposite direction at 
the same time, namely both Ei and E 2 are non-zero. In typical distribution networks, the level of voltage is the same, say 
llOV at the user level. Therefore a much more complex electric network with the same level of voltage sources can be 
equivalently transformed into the network of two voltage sources connected in the opposite direction with a series connected 
to an equivalent resistor. In the following proposition, we show that the LCL will increase as the imbalance increases between 
the two resistors of the equivalent circuit. 

Proposition 2: Suppose that Ei = E 2 = E. The LCL of Def. 1 will increase monotonically as |i?i — i? 2 | increases. 

Proof: According to Kirchhoff’s circuit Law, the current ij flows through the resistors are 


E[R2 + 2R3) 

R 1 R 2 + R 2 R 3 + R 3 R 1 

E{Ri + 2R3) 

R 1 R 2 + R 2 R 3 + R 3 R 1 
E{Ri - R 2 ) 
R 1 R 2 + R 2 R 3 + R 3 R 1 


The total system loss is 


Loss = ifRi + ilR 2 + *3f?3. 


(6) 


Since the total system loss before connecting the resistor R 3 is 


Loss' = 


{2Ef 

Ri + R2 


(7) 


we will have the LCL as follows 

Loss _ _ {Ri - R 2 Y _h 1 > 1 

Loss' 4(i?ii?2 "f R 2 R 3 "f R 3 R 1 ) 


( 8 ) 


where equality occurs if and only if Ri = i? 2 - For an unbalanced resistor network (i.e. Ri f R 2 ), rewrite the expression 
for LCL in ([^ by introducing the positive variable h = |i?i — i? 2 |- If we rewrite ([^ in terms of /i, i?i, R 3 or h, R 2 , R 3 , a 
straightforward argument using elementary calculus and counting cases shows that for any i?i, i? 2 , i ?3 > 0, i?i f R 2 , LCL 
is an increasing function of L ■ 

In decision theory, it is frequently desirable to base decisions on criteria that are known (or can be proven) to be monotonic 
in the decision variables. The previous section showed that if there is a lack of balance in the voltage source distribution 
that large current imbalances (congestion) can occur. In this section, we note that if there are imbalances in the resistances 
in our model network, then there will be a non-monotonic dependence of the losses i^Ri and i 2 i ?2 on the magnitude of 
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the resistance R 3 . 

Proposition 3: Referring to the circuit of Fig. suppose i?i < i? 2 - If Ei = E 2 , then there is a non-monotonic dependence 
of the losses i^Ri and t 2 i ?2 on the cross-link resistance R^. Specifically, there is a critical value Rf such that for R^ < R^, 
i\Ri < i 2 R 2 , while for i ?3 > R"^, i\Ri > j|i? 2 - 
Proof: We compare 

i\Ri = {E2R3 -f ii/i(i?2 + R3)) Ri/De with 
i^R^ = {E1R3 + E2{Ri + Ri/De 


where = {R1R2 + R1R3 + The loss will be larger than i|i ?2 precisely when the numerators of these 

expressions have the same magnitude relationships. Recalling Ei = E 2 = (say) E, we have 

ilRi ~ RiE^ {Rl + 4i?2i?3 + 4i?^) (9) 


and 

ilR 2 ~ R 2 E^ {Rl + AR 3 R 3 + ARl) 


( 10 ) 


Obviously, for small R 3 > 0, the expression (j^ is greater than the expression ( [TOl l, but as R 3 becomes larger, the terms 
that are quadratic in R 3 dominate—making the expression ( 101 the larger. ■ 


IV. Optimal power-flow sensitivity to small changes in network parameters 


The kinds of sensitivity illustrated in the resistive load circuits in the preceding sections may be found as well in optimal 
power flow. We revisit Example 6.16, pp. 252-254 in [5]. 

Example 1: Consider the three-node network of Fig. We examine a DC power flow model of the three bus network in 
which bus 1 and bus 2 are generators, while bus 3 is a load. The production costs of operating the generators are Cj{Pj) for 
j = 1, 2, where Pj is the nodal power injection at the j-th bus. These costs of generation are convex functions, reflecting the 
fact that as the power increases, the incremental cost rises superlinearly due to wear and tear on the machinery, decreased 
efficiency margins. The “elastic price” load at bus 3 is Pjy = Pq + Pe, where Pq is the inelastic component of the load 
and Pe is the “price elastic” component of the load. The line inductive reactances are xi 2 , 2 : 13 , T 23 . The power flow 
on line ij is given by 

H — o - 

( 11 ) 


6i — \ 


where 9j is the power phase angle at the j-hus. 
conductance matrix: 


The nodal power injections are related to the power phase angles by the 



( 12 ) 


where 



and where the line conductances 6 ^ are the negative reciprocals of the line reactances, i.e. The nodal power injections 

always sum to zero, as do the columns and rows of B, and since the power flow equations are invariant under a common 
phase shift of the Ofs it is convenient to choose a reference bus (say bus 1 ) at which we set the phase angle = 0 . 

The optimal power flow problem is to determine the nodal power injections Pi and P 2 at the generator buses and the line 
flows Pij that optimize an objective function that accounts for generation costs C'i(Pi) and C 2 {P 2 ) along with a consumer 
welfare cost that in the simplest formulation is evaluated only in terms of the price elastic load at bus 3, Cw{Pe)- The 
objective function to be minimized is written as Ci(Pi) + C 2 {P 2 ) — Cw{Pe)- Minimization is subject to constraints that 
the power flow solution does not exceed the rated capacities of the lines or buses. Thus, feasible solutions must satisfy 


0 < P, < and 0 < P„ < 


(13) 
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In a power network, congestion is said to occur if the scheduled or desired power flow exceeds the rated capacity of either 
one or more of the lines or one or more of the generator buses. 



Fig. 5. From Gomez-Exposito et al, [5]. 


Proposition 4: For the three-bus network depicted in Fig. if the power flow solution (the phases in e, @) that 
minimizes the objective function 

Ci{Pi) + C2 {P2 )-Cw{Pe) (14) 


does not result in congestion—i.e. if none of the constraints ( [T3] l holds with equality—then the optimal value of ( [l4j l is 
independent of the conductance matrix B. 

Proof: The power injections are related to the power phase angles by 


where 



(15) 


is the reduced conductance matrix, and where we have assumed without loss of generality that the power phase angle at 
bus 1 is zero. If no two of the conductances hij are zero, then is nonsingular. Using the invertible relationship (15i, we 
may rewrite the objective function as 


/(02,03)=Cl(02,03) + 4(02,03)-(?t^(02,03)- 


Solving the critical point equations df /dOj = 0 for j = 2,3 leads to a minimizing solution, and this determines the optimal 
(Pi,P 2 ) via This solution must be the same as what would have been obtained by minimizing (p^ directly. 


Remark 1: It holds much more generally that for a connected power grid, the optimal power injections from generators 
is independent of the line conductances provided there is no congestion. The proof is a direct extension of the above but is 
omitted due to space limitations. We instead examine the sensitivity of congestion to generation cost in the simple example 
at hand. 

Adopting the cost functions of [5], Cj{Pj) = PjPf and Cw{Pe) = cuPe, and recalling that Pe = Pi + P2 — Pc, the 
optimal power injections are easily seen to satisfy 

77 = ^. i = l,2. and 5 (| + £ - 2Pa) ■ 


( 16 ) 














Clearly, the cost coefficients a,/ 3 i ,/32 must be such that the power injections are within the ranges (13i—specifically the 
marginal value of consumer preference for load price elasticity must be in balance with the marginal costs of generation. 
To evaluate the line loading produced by the power injections ([T6|), we solve (15 i and use ([TT]l. This yields 


p Q ^ 12^23 _ Q ^12^12 

-^12 — 2/3i D 2^2 D 

p _ a (^12+^23)^13 I a 612^)13 

''IS — 2^1 D 2/32 D 

P a ^12^23 I a (bl2+t)13)&23 

''23 — 2I3i D 2^2 D 

where 

D — 612&13 + ^23^13 + &12^23- 


(17) 


It is expected that power grids will exhibit the same kinds if sensitivity to changes in network topology and operating 
parameters that were noted in Sections III and IV. For the generation cost values considered in Example 6.16 in [5] 
(/3i = 1, /32 = 1.675), the explicit form of the uncongested optimum power injection at generator nodes 1 and 2 favors 
power produced by the cheaper generator (generator 1), although it is never the case that the DC load flow results in zero 
power being injected at bus 2. The line loading between the less costly generator and the load (P13) turns out to be a 
monotonically increasing function of the conductance & 13 . It is interesting to note that if 613 is small enough in relation to 
623, the line loading will have P23 > 7 ^ 13 - and the load will draw more power from the line to the more expensive generator 
G 2 . 

Write the phase and line-loading relationship in matrix form; Pune = H ■ { 02 , 83 )"^, where Pum = {P 12 , P 13 , P 23 )"^ and 
H is the matrix representation specified by ([n); 



H = 


We can then express the line loadings directly in terms of the power injections by writing 


= H B- 


For small values of | 6 i 3 |, this relationship is 




( 

P12 \ 




Pl 3 



\ 

P23 

ip 

is 



/ 

P12 

\ 

( 


Pl 3 


= 

V 

P23 

) 

V 




where e = 0 ( 613 ). With the power injected by generator Gi being shifted from line (1,3) to lines (1,2) and (2,3), it is 
reasonable to expect that congestion on these lines will be sensitive to changes in Pi and generation cost parameter /3i. 
Indeed a straightforward calculation shows that there is extreme sensitivity to the cost parameter with 


P23 = 


G 

Ji 


/(|('13|), 


where G is a positive constant, and / is a smooth function of | 6 i 3 |. 

Once an uncongested optimum power flow lies outside the operating range of any component ( [T3] l, the operating limit 
of that component becomes a binding constraint in terms of which the optimal power flow problem must be resolved. 
(See [5]). Rather than pursuing constrained optimal power flow at this point, we briefly explore the pervasiveness of the 
Kirchhoff-Braess phenomena in larger networks. 


V. The Case of Large Networks 

We shall consider the effect of attaching an arbitrary two-port voltage controlled circuit (e.g. a single resistor or single 
voltage source in the simplest cases) to any two points of an existing voltage controlled circuit of arbitrary topology. It 
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will be shown that the LCL resulting from the attachment will be > 1 in all cases. We begin by recalling that a voltage 
controlled DC circuit is made up of resistors, capacitors, inductors, and voltage sources. We have the following: 

Definition 2: Points in the circuit at which two or more circuit elements are connected are called nodes. 

Fig. 1 shows a simple circuit with 5 nodes. For a DC circuit that is comprised purely of voltage sources and resistors, we 
can first choose one node to be the reference, usually called ground node. Then Kirchhoff’s Law allows for the calculation 
of the voltages at each node of the circuit, relative to the reference ground node. Once all the node voltages are known, all 
currents in the circuit can be determined easily. 



Fig. 6. A circuit with 5 nodes. 

Definition 3: For the graph of a DC circuit comprised purely of voltage sources and resistors, a fundamental node basis 
is a maximal set of nodes among which there exist no paths comprised purely of voltage source links. 

A fundamental node basis may be formed as follows. We first find the set of all nodes V where two or more circuit 
elements meet, and then create a sub-graph of the original circuit Gy = {V,Ey) that is comprised of the entire set of nodes 
V and all voltage source links. If Gy = {V,Ey) is a connected graph, any single node constitutes a fundamental node basis. 
In this case, the LCL of the circuit resulting from connecting the external two-port circuit will be equal to 1 because the 
voltage drop between any two nodes in the original circuit will be unaffected. 

If Gy = {V,Ey) is disconnected, we then can find all its connected components (with some components possibly consisting 
of single nodes). Then a fundamental node basis may be formed from the set of connected components of the given graph, 
by arbitrarily selecting one node from each connected component. Since there exists at least one path comprised purely by 
voltage source links between every node pair in a connected component, we can easily get the voltage difference between the 
fundamental node and other nodes in that connected component by computing the algebraic sum of voltage sources on their 
connecting paths. Thus once the fundamental node voltages are known, all node voltages of the circuit can be determined 
easily without solving the Kirchhoff’s equations. 

Proposition 5: If a DC circuit is comprised purely of voltage sources, resistors, capacitors and/or inductors, then the 
steady state LCL resulting from ADDING a resistance link must be > 1 no matter what topological structure the original 
circuit has. 

Proof: We start with the simplification of the circuit. If a DC circuit comprised purely of voltage sources, resistors, 
capacitors and/or inductors is in its steady state, then we can remove all loops with one or more capacitors (since capacitors 
act identically as open circuits in DC steady state) and replace all inductors with short circuits (since inductors act identically 
as short circuits in DC steady state) while keeping the total loss of the circuit unchanged. 

We then can find all nodes where two or more circuit elements meet (only resistors or voltage sources are left after the 
above simplification), and mark them as V = {1,..., /V}, where N > 2 is the total number of nodes. Then each link between 
any pair of adjacent nodes {i,j} {i,j = 1, and i j) will only have either one voltage source or one resistor. 

We next create the voltage source sub-graph of the original circuit Gy = {V,Ey) which is comprised of the whole set 
of nodes V = iV} and all voltage source links. The case that the voltage source sub-graph Gy = {V,Ey) is a 
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connected graph has been seen to be trivial and thus is ignored here. Suppose Gy = {V,Ey) is disconnected, we then can 
find all its connected components and assume that they are Gi = {Vi, Ei),G m = (Ym, Em) (M is the total number of 
connected components). A fundamental node basis then can be formed by arbitrarily selecting one node from each connected 
component, and we denote this basis by Vp = {r'l, ...,um} (1 < vi, ...jVm < E). We denote the voltage at node Vi by 
(i = 1, ...,M), and we choose vm as the reference ground node, i.e. ey,^, = 0. 

It is easy to prove that the pair of endpoints of any resistor link in the original circuit belong to either one connected 
component or two connected components defined above. For those resistor links whose pair of endpoints belong to one 
connected component, their power loss will be unchanged no matter what kind of new link is introduced as the voltage 
drops between the pairs of endpoints are fixed. Thus we can use a scalar constant, say Pi, to denote the total loss of such 
resistor links. 

Next we compute the loss of the resistor links whose pair of endpoints belong to two connected components. We denote 
the number of resistor links between the pair of connected components {Gi,Gj} {i,j = 1, ...,M,and i ^ j) by Lij, the 
current flowing away from the i-th connected component towards the j-th connected component on the /c-th resistor link 
between {Gi,Gj} {k = by the resistor on the fc-th link between {Gi,Gj} by Rij^k 

iRi,j,k=Rj,i,k > 0). Then the total loss of such resistor links is given by 

M-l M Lij 

(18) 

i—l j— 2+1 k —1 


where hj^k can be expressed as 


k,j,k 


ep — Cy ■ — ep 


R 


i,j,k 


(19) 


with ep^ ^ denoting the algebraic sum of voltage sources on the path connecting the fundamental node Vi and the node 
to which the resistor Rij,k is attached in the connected component Gi. (ep„ ^ = 0 if one endpoint of the resistor Ri,j,k 
is directly connected to the fundamental node Vi.), Similarly, ep is the sum of voltages along the path connecting the 
other endpoint of the resistor to the fundamental node Vj. (ep„ ^ =0 if the other endpoint of the resistor Ri,j,k is directly 
connected to the fundamental node vj.) 

A potential function whose physical meaning is the total loss of all resistors in the original circuit can be created based 
on the variables in the fundamental node set. 

This is given explicitly by 


P = Pi+P 2 

M-l M Lij , . _p 

i=l j=i+l k=l 


( 20 ) 


where {et,^,..., ey,^,} are the node voltages. We assume = 0, and the value of other node voltages before and after adding 
a new link are (e„j,..., and ..., J, respectively. We shall show that the potential function reaches its 

minimum loss level before adding a new link, i.e. ,..., is always < P{e ^_^, ..., 

We first show that the potential function is a strict convex function. Suppose that 


„ (e„. + ep„^_^ - - ep„^_ 

R,j,k — P 


Pi,j,k is obviously a convex function of the node voltages. Since we assume = 0 and suppose 

{cy. + ep - ep„,^ 


P 


i.M.k — 


Ri,M,k 

then Pi^M.k is also strictly convex provided Ri^mY finite- Since 




M-l M Li^j 

p = p^ + Jl E E^^^a 

2 = 1 j = 2+l k—1 


( 21 ) 


( 22 ) 


P is also a strictly convex function, taking its global minimum where all its partial derivatives are zero. 


(23) 
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We show that all partial derivatives of potential function are zero when , et,^,J = (e„j,J. Since 


dP 


i —1 

j=i k=i 
M Lij 

EE 

k—1 


ep„ . _ — e„ — ep 


’j,fc 


^vi + , — G.V- — ep^_ 


R 


i,j,k 


(24) 


the partial derivative of P in the direction e„. is exactly double the algebraic sum of all currents flowing on the original 
resistors that meet at the *-th connected component. Because the node voltages 

must satisfy Kirchhoff’s Current Law, the partial derivatives in Equ. |^must be zero for j = 1,... ,M — 1. The potential 
function thus reaches its global minimum loss level under normal operating conditions. We consider what happens if a new 
link (two-port circuit) is attached to any pair of nodes. If both endpoints of the newly introduced link are added to the 
same connected component of Gy(V,Ey), defined above, then the total loss of the original circuit will remain the same as 
(e„^,..., e„„_ J J. However, if the two endpoints of the new link are added to two different connected 

components, say Gi and Gj, and if we denote the current flowing on the new link from Gi to Gj by I"®’", then we will 
have 


dP 

de'. 


2—1 Lij 

j^l k^l 
M 

E 

2 + 1 k—1 

= 2 X /"+. 


E' 


^Vi + - ep„ 

Ri,j,k 

+ ep„ . , — e — ep 




(25) 


Ri,j,k 


After adding the new link, the algebraic sum of all currents flowing on the original resistors that meet at the i-th connected 
component is determined by the current 7”®“ which is not necessarily zero as there may be current import and export to 
the newly added link from the original system. Therefore the original loss equilibrium is perturbed away from its minimum 
level, and thus the system loss of the original system will have increased. This proves the proposition. ■ 

Remark 2: The multi-node connected components of Gy(V,Ey) may have arbitrarily complex topologies-including tree 
and loop components. If a fundamental node basis that differs from the one chosen to define the loss in p0| ) is chosen, the 
form of the loss function (201 will differ accordingly. The critical point determined by setting the partial derivatives in (24 1 
equal to zero will minimize the new expression for loss. It follows from the Kirchhoff circuit laws, that the minimizing 
values in both representations are the same—as we would expect. 

Remark 3: A more general version of Proposition]^ can be established. Indeed, the proof as given applies to the connection 
any two-port voltage controlled circuit to an existing voltage controlled circuit in steady state. Although the proof becomes 
more involved, a similar result holds for the addition of an n-port external circuit. 


VI. Conclusions 

It has long been recognized that contingencies like the loss of a major power line can pose significant threats to the secure 
operation of the power grid. The focus of this paper has been on the way that seemingly small changes can have large 
effects. We have presented examples of simple circuits and networks that display the kinds of sensitivity to small changes 
in operating parameters that will need to be better understood as smart microgrids become an increasingly important part 
of power distribution networks. Secure operation of these microgrids will require the real-time coordinated control of 
increasing numbers of small-scale generation resources and consumer-provided demand response assets while respecting the 
safe operating ranges of all lines and equipment. Taking inspiration from concepts of congestion in traffic networks, we have 
studied what we call Kirchhoff-Braess phenomena—the apparent worsening of congestion due to the addition of a lightly 
loaded line connecting points in the network. We have defined a power network analogue of the price of anarchy that we 
call the Loss Cost of the Link (LCL). This is simply the ratio of losses after and before the asset (e.g. a line) has been 
added. It has been shown that this ratio is frequently greater than one, and this has been studied in detail for simple small 
networks and in considerable generality for voltage controlled networks in Section VI. Space does not permit treatment of 
current-controlled circuits, but a corresponding theory characterizes situations in which network losses increase and other 
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situations in which there are decreases in network losses when a new circuit link is added. An important goal of future 
research is an understanding of general classes of optimal power flows in which the cost sensitivity of network congestion 
as treated in Section V is extended to realistic transmission and distribution networks. Models that capture degradation of 
voltage, frequency and other important physical parameters are being developed. 
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